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1.  INTRODUCTION 


Sprays  and  spray  processes  have  been  studied  extensively  due  to  their  many 
applications  (Faeth  1987,  1990).  Nevertheless,  sprays  are  complex  multiphase  turbulent 
flows  and  fundamental  understanding  of  their  properties  is  not  well  developed  — 
particularly  those  processes  relating  to  die  near-injector  dense  portion  of  the  flow.  This  is 
a  serious  deficiency  because  the  dense-spray  region  involves  the  breakup  of  the  liquid  to 
form  a  dispersed  phase,  which  is  crucial  to  the  mixing  properties  of  sprays,  and  also 
develops  the  initial  conditions  required  to  analyze  the  structure  of  the  better  understood 
dilute  portion  of  sprays.  Thus,  the  objective  of  the  present  investigation  was  to  study 
three  aspects  of  dense  sprays  that  involve  drop/gas  interactions,  namely:  (1)  secondary 
drop  breakup,  which  is  an  intrinsic  outcome  of  primary  breakup  of  the  liquid,  and  is  the 
most  significant  rate  process  of  dense  sprays;  (2)  turbulence  generation  by  dispersed 
phases,  which  is  the  most  significant  source  of  turbulence  production  within  dense 
sprays;  and  (3)  the  structure  of  sphere  wakes  at  moderate  Reynolds  numbers,  which  is  a 
fundamental  property  needed  to  understand  the  mechanism  of  turbulence  generation.  The 
research  has  relevance  to  air-breathing  propulsion  systems,  liquid  rocket  engines  and 
internal  combustion  engines,  among  others. 

The  following  description  of  the  research  is  relatively  brief.  Additional  details 
may  be  found  in  articles,  papers  and  theses  resulting  from  the  investigation  that  are 
summarized  in  Table  1.  This  table  also  provides  a  summary  of  participants  in  the 
investigation  and  oral  presentations  of  portions  of  the  research  results.  Finally,  for 
convenience,  several  articles  resulting  from  the  research  are  reproduced  in  appendices, 
including:  Faeth  (1990),  Hsiang  and  Faeth  (1992),  Parthasarathy  and  Faeth  (1990a,b), 
Mizukami  et  al.  (1992)  and  Wu  and  Faeth  (1992). 

Table  1.  Summary  of  Investigation 


Articles.  Papers  and  Theses: 

Faeth,  G.  M.  (1990)  Structure  and  atomization  properties  of  dense  turbulent  sprays. 
Twentv-third  Symposium  (International)  on  Combustion.  The  Combustion 
Institute,  Pittsburgh,  1345-1352. 

Faeth,  G.  M.  and  Parthasarathy,  R.  N.  (1989)  Homogeneous  particle- laden  flows:  2. 
turbulent  dispersion.  Bull.  Amer.  Phvs.  Soc.  34, 2312  (abstract  only). 

Faeth,  G.  M.  and  Wu,  J.  S.  (1991)  Structure  of  the  wakes  of  spheres  at  moderate 
Reynolds  numbers.  Bull.  Amer.  Phvs.  Soc.  36,  2625  (abstract  only). 

Faeth,  G.  M.,  Mizukami,  M.  and  Parthasarathy,  R.  N.  (1990)  Particle-generated 
turbulence  in  homogeneous  dilute  dispersed  flows.  Bull.  Amer.  Phvs.  Soc.  35, 
2286  (abstract  only). 

Hsiang,  L.-P.  and  Faeth,  G.  M.  (1992)  Near-limit  drop  deformation  and  secondary 
breakup.  Int.  J.  Multiphase  Flow  1 8, 635-652. 

Hsiang,  L.-P.  and  Faeth,  G.  M.  (1993)  Deformation  and  secondary  breakup  of  drops. 
AIAA  Paper  No.  93-0814. 
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Mizukami,  M.,  Parthasarathy,  R.  N.  and  Faeth,  G.  M.  (1992)  Particle-generated 
turbulence  in  homogeneous  dilute  dispersed  flows.  Int.  J.  Multiphase  Flow 
18,397-412. 

Parthasarathy,  R.  N.  (1989)  Homogeneous  dilute  turbulent  particle-laden  water  flows. 
Ph.D.  thesis,  The  University  of  Michigan,  Ann  Arbor,  Michigan. 

Parthasarathy,  R.  N.  and  Faeth,  G.  M.  (1989)  Homogeneous  particle-laden  flows:  1. 
continuous-phase  turbulence  properties.  Bull.  Amer.  Phvs.  Soc.  34,  2312  (abstract 
only). 

Parthasarathy,  R.  N.  and  Faeth,  G.  M.  (1990a)  Turbulence  modulation  in  homogeneous 
dilute  particle-laden  flows.  J.  Fluid  Mech.  220, 485-514. 

Parthasarathy,  R.  N.  and  Faeth,  G.  M.  (1990b)  Turbulent  dispersion  of  particles  in  self¬ 
generated  homogeneous  turbulence.  J.  Fluid  Mech.  220, 515-537. 

Ruff,  G.  L.  (1990)  Structure  and  mixing  properties  of  the  near-injector  region  of 
nonevaporating  pressure-atomized  sprays.  Ph.D.  thesis,  The  University  of 
Michigan,  Ann  Arbor,  Michigan. 

Ruff,  G.  A.,  Bernal,  L.  P.  and  Faeth,  G.  M.  (1990)  High  speed  in-line 
holocinematography  for  dispersed-phase  dynamics.  J.  AddI  Optics  29,  4544- 
4546. 

Ruff,  G.  A.,  Bernal,  L.  P.  and  Faeth,  G.  M.  (1991)  Structure  of  the  near-injector  region  of 
non-evaporating  pressure-atomized  sprays.  J.  Prop.  Power  7. 22 1  -230. 

Ruff,  G.  A.,  Wu,  P.-K.,  Bernal,  L.  P.  and  Faeth,  G.  M.  (1992)  Continuous-and  dispersed- 
phase  structure  of  dense  nonevaporating  pressure-atomized  sprays.  J.  Pron.  Power 
8, 280-289.  . 

Wu,  J.  S.  and  Faeth,  G.  M.  (1992a)  Sphere  wakes  in  still  surroundings  at  intermediate 
Reynolds  numbers.  AIAA  J..  in  press. 

Wu,  J.  S.  and  Faeth,  G.  M.  (1992b)  Sphere  wakes  in  turbulent  surroundings  at 
intermediate  Reynolds  numbers.  AIAA  J..  submitted. 

Participants: 

Faeth,  G.  M.,  Principal  Investigator,  Professor,  The  University  of  Michigan 

Hsiang,  L.-P.,  Graduate  Student  Research  Assistant,  Doctoral  Candidate,  The  University 
of  Michigan. 

Mizukami,  M.,  Graduate  Student  Research  Assistant,  The  University  of  Michigan; 
currently.  Research  Engineer,  Propulsion  Systems  Division,  NASA  Lewis 
Research  Center. 

Parthasarathy,  R.  N.,  Graduate  Student  Research  Assistant,  Doctoral  Candidate,  The 
University  of  Michigan;  currently  Research  Scientist,  Institute  of  Hydraulic 
Research,  University  of  Iowa. 
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Ruff,  G.  A.,  Graduate  Student  Research  Assistant,  Doctoral  Candidate,  The  University  of 
Michigan;  currently,  Assistant  Professor,  Mechanical  Engineering,  Drexel 
University. 

Wu,  J.  S.,  Graduate  Student  Research  Assistant,  Doctoral  Candidate,  The  University  of 
Michigan. 

Oral  Presentations: 

Faeth,  G.  M.,  "Separated-flow  and  breakup  phenomena  in  dense  sprays,"  ASME 
Workshop  on  Issues  in  Spray  Systems,  Philadelphia,  PA,  August  1989. 

Faeth,  G.  M.,  "Structure  of  dense  pressure-atomized  sprays, "Mechanical  Engineering 
Seminar,  Camegie-Mellon  University,  Pittsburgh,  PA,  September  1989. 

Faeth,  G.  M.,  "Homogeneous  particle-laden  flows.  2.  Turbulent  dispersion,"  42nd  Annual 
Meeting,  American  Physical  Society,  Division  of  Fluid  Dynamics,  Palo  Alto,  CA, 
November  1989. 

Faeth,  G.  M.,  "Progress  concerning  dense  spray  processes,"  International  Atomization 
and  Spray  Systems  Conference,  Hartford,  CT,  May  1990. 

Faeth,  G.  M.,  "Atomization  and  mixing  in  dense  turbulent  sprays,"  Twenty-Third 
Symposium  (International)  on  Combustion,  Orleans,  France,  July  1990. 

Faeth,  G.  M.,  'Turbulence  generation  in  homogeneous  dilute  particle/liquid  or  gas  flows," 
Mechanical  Engineering  Seminar,  Wayne  State  University,  Detroit,  MI, 
September  1990. 

Faeth,  G.  M.,  "Particle-generated  turbulence  in  homogeneous  dilute  dispersed  flows," 
43rd  Annual  Meeting,  American  Physical  Society,  Division  of  Fluid  Dynamics, 
Ithaca,  NY,  November  1990. 

Faeth,  G.  M.,  "Structure  and  breakup  processes  of  dense  sprays,"  Mechanical  and 
Aerospace  Engineering  Seminar,  Princeton,  NJ,  December  1990. 

Faeth,  G.  M.,  "Optical  diagnostics  for  combustion  studies,"  NSF  Workshop  on  Particulate 
Two-Phase  Flow  Visualization,  Washington,  DC,  March  1991. 

Faeth,  G.  M.,  "Structure,  breakup  and  turbulence  interactions  in  sprays,"  NSF  Workshop 
on  Atomization  and  Sprays  2000,  Gaithersburg,  MD,  July  1991. 

Faeth,  G.  M.,  "Structure  of  the  wakes  of  spheres  at  moderate  Reynolds  numbers,"  44th 
Annual  Meeting,  American  Physical  Society,  Division  of  Fluid  Dynamics, 
Tempe,  AZ,  November  1991. 

Faeth,  G.  M.,  "Secondary  drop  breakup  in  the  deformation  regime,"  AIAA  30th 
Aerospace  Sciences  Meeting,  Reno,  NV,  January  1992. 

Hsiang,  L.-P.,  "Near-limit  drop  breakup  from  step  vel  xtity  changes,"  International 
Conference  on  Liquid  Atomization  and  Spray  Systems,  Gaithersburg,  MD,  July 
1991. 
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Hsiang,  L.-P.,  "Secondary  drop  breakup,"  Gas  Dynamics  Seminar,  University  of 
Michigan,  Ann  Arbor,  MI,  October  1991. 

Parthasarathy,  R.  N.,  "Homogeneous  particle-laden  flows:  1.  Continuous-phase 
turbulence  properties,"  42nd  Annual  Meeting,  American  Physical  Society, 
Division  of  Fluid  Dynamics,  Palo  Alto,  CA,  November  1989. 

Wu,  J.-S.,  "Structure  of  sphere  wakes  at  moderate  Reynolds  numbers,"  Gas  Dynamics 
Seminar,  University  of  Michigan,  Ann  Arbor,  MI,  February  1992. 


The  following  report  considers  secondary  breakup,  turbulence  generation,  and 
sphere  wakes,  in  turn.  Each  section  is  written  so  that  it  stands  alone;  therefore,  readers 
can  skip  to  sections  of  interest. 


2.  SECONDARY  BREAKUP 


2.1  introduction 

Secondary  breakup  of  drops  is  an  important  fundamental  multiphase  flow 
phenomenon  with  applications  to  propellant  combustion  and  numerous  other  natural  and 
engineering  processes.  In  particular,  recent  studies  of  dense  sprays  confirm  the 
conventional  view  of  liquid  atomization  that  primary  breakup  at  a  liquid  surface  yields 
drops  that  are  intrinsically  unstable  to  secondary  breakup;  and  that  secondary  breakup 
tends  to  control  mixing  rates  in  dense  sprays  much  like  drop  vaporization  tends  to  control 
mixing  rates  in  dilute  sprays,  see  Ruff  et  al.  (1989a, b,  1991,  1992)  and  Faeth  (1987, 
1990)  and  references  cited  therein.  Motivated  by  these  observations,  the  objectives  of  this 
phase  of  the  investigation  were  to  study  drop  deformation  and  breakup  for  well-defined 
shock  wave  disturbances. 

Reviews  of  past  work  on  secondary  breakup  are  reported  by  Hinze  (1955), 
Krzeczkowski  (1980),  Faeth  (1987,  1990)  and  Hsiang  and  Faeth  (1992).  The  reviews 
indicate  that  the  regimes  of  breakup  have  been  defined  reasonably  well  but  analogous 
deformation  regimes  are  unknown,  that  effects  of  liquid  viscosity  (high  Oh)  on  breakup 
times  have  not  been  examined  and  that  information  concerning  the  outcome  of  breakup  is 
not  available.  Thus,  these  issues  were  emphasized  during  the  present  theoretical  and 
experimental  investigation.  The  present  description  of  the  research  is  brief,  see  Hsiang 
and  Faeth  (1992, 1993)  for  more  details. 

2-2  Experimental  Methods* 

Apparatus.  A  shock  tube  with  the  driven  section  open  to  the  atmosphere  was  used 
to  generate  shock- wave  disturbances,  see  Fig.  1.  The  driven  section  was  rectangular  (38  x 
64  mm)  with  a  6.7  m  length  to  provide  17-24  ms  of  uniform  flow  behind  the  shock  wave 
for  test  purposes.  A  vibrating  capillary  tube  generator  provided  a  stream  of  drops  across 
the  shock  tube,  with  an  electrostatic  separation  system  used  to  vary  spacing  between 
drops  in  order  to  avoid  drop/drop  interactions.  The  interactions  between  the  drop  stream 
and  the  shock  wave  flow  could  be  observed  through  windows  in  the  walls  of  the  driven 
section  (25  mm  high  x  305  mm  long). 

Instrumentation.  Drops  were  observed  in  two  ways:  pulsed  shadowgraph 
photographs  and  motion  pictures  to  observe  overall  breakup  dynamics,  and  single-and 
double-pulse  holography  to  observe  the  outcome  of  breakup.  The  details  of  these 
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igure  1.  Sketch  of  the  shock  tube  apparatus 


instruments  and  the  resulting  uncertainties  of  the  measurements  are  discussed  by  Hsiang 
and  Faeth  (1992, 1993) 

2.3  Results  and  Discussion 


The  presentation  of  results  will  begin  with  definition  of  deformation  and  breakup 
regime  transitions  to  help  organize  the  rest  of  the  findings.  These  transitions  are  plotted  in 
Fig.  2  as  functions  of  We  and  Oh,  following  earlier  work  (Hinze  1955;  Krzeczkowski 
1980).  Three  breakup  regimes  are  delineated;  bag  breakup,  where  the  center  of  the  drop 
deforms  as  a  bag  in  the  downstream  direction;  shear  breakup,  where  drops  are  stripped 
from  the  periphery  of  the  original  drop;  and  multimode  breakup,  which  is  a  complex 
combination  of  the  first  two  modes.  In  view  of  the  somewhat  arbitrary  definition  of  these 
regimes,  the  agreement  of  present  results  with  those  of  Hinze  (1955)  and  Krzeczkowski 
(1980)  is  quite  satisfying.  The  deformation  regimes  illustrated  in  Fig.  2  have  not  been 
reported  before,  and  involve  non  oscillatory  and  oscillatory  deformation  with  the  latter 
regime  being  relatively  limited. 

The  most  striking  feature  deformation  and  breakup  regime  map  of  Fig.  2  is  that 
progressively  higher  We  numbers  are  needed  for  the  various  transitions  as  Oh  increases. 
In  fact,  it  appears  that  breakup  might  no  longer  be  observed  for  Oh  somewhat  greater  than 
4,  which  was  the  highest  value  observed  during  the  present  investigation,  with 
deformation  also  probably  disappearing  for  somewhat  larger  values  of  Oh.  This  large  Oh 
regime  is  encountered  during  spray  combustion  at  high  pressures  where  the  surface 
tension  becomes  small  as  the  drop  approaches  thermodynamic  critical  conditions.  Thus, 
present  findings  suggest  that  drops  will  not  shatter  into  small  drops  near  the 
thermodynamic  critical  point,  as  often  thought  (Faeth  1990);  instead,  they  will  only 
deform,  or  simply  remain  spherical  at  sufficiently  large  Oh.  This  behavior  has  significant 
implications  concerning  spray  combustion  processes  at  elevated  pressures,  however, 
resolving  the  issue  will  require  considering  liquid/gas  density  ratios  much  nearer  to  unity 
than  the  present  test  conditions. 

Drop  breakup  time  is  a  key  parameter  needed  to  understand  effects  of  liquid 
viscosity  through  Oh  on  breakup  processes,  e.g.,  as  breakup  and  velocity  relaxation  times 
approach  one  another,  the  propensity  for  drop  breakup  decreases  due  to  reduced  relative 
velocities.  Present  measurements  of  breakup  times,  along  with  earlier  results  for  shock 
wave  disturbances  due  to  Engel  (1958),  Simpkins  and  Bales  (1972),  Ranger  and  Nicholls 
(1969)  and  Reinecke  and  coworkers  (1969,  1970),  are  plotted  as  a  function  of  We  in  Fig. 
3.  The  breakup  times  are  normalized  by  the  characteristic  shear  breakup  time  of  Ranger 
and  Nicholls  (1969): 


t*  =  do(piypG)1/2/uo  (1) 

Except  for  present  results,  which  are  grouped  according  to  Oh,  the  measurements  are  for 
OhcO.l  where  effects  of  liquid  viscosity  are  small.  A  remarkable  feature  of  the  results  is 
that  tb/t*  varies  very  little  even  though  We  varies  from  10  to  106  for  Oh  <  0.1,  yielding 
tb/t*  =  5.  However,  tb  progressively  increases  as  Oh  increases,  yielding  the  following 
empirical  relationship  over  the  present  test  range 

tb/t*  =  5/(  1  -Oh/7);  We  <  103  (2) 

This  result  is  consistent  with  the  observation  of  a  limited  breakup  region  as  Oh  increases 
but  is  based  on  relatively  few  data  with  Oh  <  3.5. 
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Figure  3.  Drop  breakup  times  as  a  function  of  We  and  Oh 


Similar  to  past  work  on  the  structure  of  dense  sprays  and  processes  of  primary 
breakup  (Ruff  et  al.  1992;  Wu  et  al.  1992c),  drop  size  distributions  after  secondary 
breakup  generally  satisfied  Simmons'  (1977)  universal  root  normal  distribution  with 
MMD/SMD  =  1.2.  (Hsiang  and  Faeth  1992).  Shear  breakup  was  an  exception  when  the 
drop-forming  drop  was  considered,  however,  this  deficiency  could  be  eliminated  by 
removing  the  drop-forming  drop  from  the  distribution  and  treating  it  separately  (Hsiang 
and  Faeth  1993).  Thus,  the  entire  drop  size  distribution  can  be  represented  by  the  SMD 
alone,  which  is  the  practice  that  will  be  followed  subsequently. 

A  correlation  for  the  sizes  after  secondary  breakup  was  developed  based  on 
phenomenological  analysis  of  shear  breakup  that  served  equally  well  for  the  bag  and 
multimode  breakup  processes  (Hsiang  and  Faeth  1992).  Present  measurements  of  drop 
sizes  after  secondary  breakup  are  illustrated  according  to  this  correlation  in  Fig.  4.  These 
results  are  for  Oh  <  0.1  and  We  <  103,  and  yield  the  following  fit: 

pG  SMD  Uo2/a  =  6.2  (piypG)1/4  [PLftpL  do  Uq)]1/2  We  (3) 

with  the  overall  correlation  coefficient  of  the  fit  being  0.91.  It  should  be  noted,  however, 
that  Pl/pg  does  not  vary  greatly  over  the  present  list  range  and  additional  measurements 
are  needed  to  explore  density  ratio  effects.  It  probably  is  fortuitous,  and  certainly 
surprising,  that  a  single  correlation  can  express  the  SMD  after  bag,  multimode  and  shear 
breakup.  However,  this  behavior  is  consistent  with  the  observation  that  their  breakup 
times  correlate  in  the  same  way,  see  Fig.  3. 

The  breakup  regimes  of  Fig.  2  are  marked  on  Fig.  4,  assuming  that  the  velocities 
of  the  drops  after  secondary  breakup  remain  at  u<>.  If  this  were  the  case,  drops  at  initially 
large  We  would  remain  unstable  to  breakup  after  the  secondary  breakup  process  had 
occurred,  which  obviously  is  not  plausible.  Recent  work  has  resolved  this  difficulty  by 
finding  drop  velocity  distributions  after  secondary  breakup  and  successfully  correlating 
these  results  using  phenomenological  theory  (Hsiang  and  Faeth  1993).  This  shows  that 
variations  of  drop  velocities  during  secondary  breakup  are  sufficient  to  reduce  We 
numbers  to  the  stable  regime.  An  exception  is  the  drop  forming  drop  during  high  We 
shear  breakup,  however,  in  this  case,  the  drop  forming  drop  no  longer  is  subject  to  shock 
loading  so  that  it  becomes  stable  even  at  a  large  We.  These  results  highlight  the  fact  that 
We,  time  and  the  relative  acceleration  rates  of  drops  all  affect  drop  breakup,  with  shock 
loading  only  being  a  convenient  limiting  process,  see  Hsiang  and  Faeth  (1993)  for 
additional  discussion  of  this  behavior. 

2.4  Conclusions 

Drop  deformation  and  secondary  breakup  after  shock  wave  initiated  disturbances 
were  studied  for  a  variety  of  liquids  in  air  at  normal  temperature  and  pressure.  Major 
conclusions  of  the  study  are  as  follows: 

1.  The  We  for  the  onset  of  deformation  and  breakup  regimes  increase  for  increasing 
Oh,  with  no  breakup  observed  over  the  present  test  range  for  Oh  >  4.  This 
suggests  that  drops  do  not  shatter  as  they  approach  their  thermodynamic  critical 
point,  as  generally  thought  (Faeth  1990);  instead,  they  simply  deform  or  may  even 
remain  as  spheres  for  sufficiently  large  Oh. 

2.  Unified  scaling  of  breakup  times  was  observed,  showing  relatively  weak  effects  of 
We,  but  increasing  breakup  times  at  Oh  increased. 
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Figure  4.  Correlation  of  SMD  after  secondary  breakup 


10 


3.  Drop  size  distributions  after  secondary  breakup  satisfied  Simmons’  (1977) 
universal  root  normal  distribution  with  MMD/SMD  =  1.2.  The  SMD  after 
secondary  breakup  could  be  correlated  based  on  phenomenological  analysis  of 
shear  breakup  for  Oh  <  0. 1  and  all  three  breakup  regimes. 

4.  It  also  was  possible  to  correlate  drop  velocity  distributions  after  secondary 
breakup  based  on  phenomenological  analysis  of  the  shear  breakup  regime.  These 
results  show  that  most  drops  are  below  stability  limits  for  shock  disturbances  after 
secondary  breakup  while  the  drop  forming  drop  during  shear  breakup  is  stable  for 
its  accelerative  environment  when  drop  stripping  ends. 

Most  of  present  test  results  are  for  Oh  <  0.1  and  Pl/pg  >  500.  Thus,  current  work 
is  emphasizing  test  conditions  at  larger  Oh  and  smaller  Pl/Pg  »  that  are  more 
representative  of  combusting  sprays  at  the  elevated  pressures  of  contemporary  propulsion 
systems. 


3.  TURBULENCE  GENERATION 


3.1  Introduction 

The  objective  of  this  portion  of  the  investigation  was  to  study  turbulence 
generation  by  the  motion  of  dispersed  phases  in  multiphase  flows.  Experimental 
conditions  involved  uniform  fluxes  of  particles  falling  in  nearly  stagnant  (in  the  mean) 
air,  in  order  to  supplement  earlier  findings  for  similar  particle/water/flows  (Parthasarathy 
and  Faeth  1990a,b).  A  simplified  stochastic  model  of  the  flow  due  to  Parthasarathy  and 
Faeth  (1990a)  was  used  to  interpret  and  correlate  both  sets  of  measurements. 

Reviews  of  past  work  on  turbulence  generation  are  reported  by  Parthasarathy  and 
Faeth  (1990a,b)  and  Mizukami  et  al.  (1992).  Based  on  this  status,  the  objective  of  the 
present  investigation  was  to  extend  the  particle/water  study  of  Parthasarathy  and  Faeth 
(1990a)  to  particle/air  flows.  The  main  motivation  for  this  step  is  that  rates  of  dissipation 
of  particle  energy  in  air  are  orders  of  magnitude  larger  than  for  particles  in  water  at 
similar  conditions  so  that  effects  of  this  important  parameter  can  be  resolved.  The 
following  description  of  the  research  is  brief,  additional  details  can  be  found  in 
Parthasarathy  and  Faeth  (1990a,  b)  and  Mizukami  et  al.  (1992). 

3.2  Experimental  Methods 

Apparatus.  A  sketch  of  the  particle/air  flow  apparatus  appears  in  Fig.  5.  The  tests 
were  conducted  using  monodisperse  spherical  glass  particles  having  nominal  diameters  of 
0.5,  1.0  and  2.0  mm.  The  particles  were  dispersed  by  falling  through  an  array  of  screens 
and  then  passed  through  a  windowed  test  chamber  (410  x  535  x  910  mm)  where 
measurements  were  made.  The  particles  collected  with  little  rebound  at  the  bottom  of  the 
chamber  and  were  removed  periodically. 

Instrumentation.  Measurements  involved  particle  number  fluxes,  particle 
velocities  and  continuous-phase  flow  properties.  Particle  number  fluxes  were  measured 
by  collecting  particles  in  containers  at  the  bottom  of  the  test  chamber  for  timed  intervals: 
they  were  uniform  within  experimental  uncertainties  except  for  the  region  near  the  walls 
of  the  test  chamber.  Particle  velocities  were  measured  by  particle  tracking,  illuminating 
the  central  portion  of  the  tank  with  a  stroboscopic  light  source  and  recording  the  images 
with  an  open  camera  shutter.  Continuous-phase  mean  and  fluctuating  velocities  were 
measured  with  a  two-point  phase  discriminating  laser  velocimeter  (LV),  identical  to 
Parthasarathy  and  Faeth  (1990a,b).  Assessment  of  flow  properties  with  the  LV  showed 
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that  they  were  uniform  within  the  central  portion  of  the  test  chamber  (±  120  mm  in  the 
cross-stream  direction  and  ±180  mm  in  the  stream  wise  direction). 

3.3  Theoretical  Methods 

The  simplified  analysis  of  Parthasarathy  and  Faeth  (1990a)  was  used  to  help 
interpret  the  measurements.  The  major  assumptions  of  the  analysis  are  as  follows:  the 
flows  are  statistically  stationary  with  uniform  particle  fluxes  and  constant  continuous- 
phase  properties;  particle  arrival  times  are  independent  of  other  particle  arrival  times  so 
that  they  satisfy  Poisson  statistics;  the  flows  are  infinite  in  extent,  the  flows  are  dilute  so 
that  the  probability  of  a  test  point  being  within  a  particle  is  small;  the  contribution  of  flow 
properties  immediately  around  particles  is  ignored  for  the  same  reason;  and  because  the 
equations  of  motion  are  linear  for  asymptotic  wakes,  flow  properties  are  taken  to  be  the 
result  of  a  linear  superposition  of  particle  flow  fields  (wakes)  that  have  reached  a 
particular  position. 

Summing  flow  properties  under  these  assumptions  is  an  extension  of  methods 
used  to  analyze  random  noise  (Rice  1954).  Taking  the  point  of  observation  be  the  origin 
of  a  cylindrical  coordinate  system  and  taking  the  effect  of  a  particle  to  have  mean  and 
turbulent  contributions  of  g  (r,  <j>,  t)  and  g'  (r,  t),  the  resulting  mean-squared  fluctuation 

about  the  average  (which  is  zero  for  present  flows)  becomes  (Parthasarathy  and  Faeth 
1990a): 


G'2  =  n"  f~  dt  ft*  d<|>  ft  [g2  (r,  4>,  0  +  g’2  (r,  <j>,  t)]  r  dr.  (4) 

Equation  (4)  is  for  a  monodisperse  particle  or  drop  flow.  Under  the  assumption  of  linear 
superposition,  however,  if  <|>  (d)  is  a  generic  property  for  a  particular  diameter,  then  the 
mean  value  of  <(>  becomes: 


$  =  Jo  <t>(d)PDF(d)dd  (5) 

where  PDF  (d)  is  the  probability  density  function  of  particle  diameter. 

The  main  difficulty  in  applying  Eqs.  (4)  and  (5)  is  that  little  is  known  about  the 
properties  of  particle  wakes  at  modest  Reynolds  numbers  (<  100)  in  a  turbulent 
environment,  as  discussed  later.  Effects  of  wake  turbulence  cannot  be  ruled  out  because 
wake  Reynolds  numbers  are  greater  than  unity  for  the  range  of  interest;  therefore, 
integration  to  find  flow  properties  used  mean  turbulent  wake  properties  from  existing 
results  for  asymptotic  turbulent  wakes  (Uberoi  and  Freymuth  1970;  Tennekes  and 
Lumley  1972).  A  second  problem  is  that  the  integrals  of  Eq.  (4)  do  not  converge  as  t  — » 
oo.  Pending  resolution  of  this  difficulty,  integrations  were  terminated  at  z/d  =  175,  which 
provides  a  reasonable  fit  of  the  measurements  and  yields  mean  wake  velocities  that  are 
small  in  comparison  to  velocity  fluctuation  levels  which  would  be  expected  to  be  lost  in 
the  background  turbulent  field. 

3.4  Results  and  Discussion 

Particle-generated  turbulence  fields  are  unusual  because  it  is  rather  simple  to  find 
the  local  rate  of  dissipation  of  turbulence  kinetic  energy,  which  must  be  equal  to  the  local 
loss  of  particle  mechanical  energy,  e.g. 

e  =  rc  n"  d2  CD  U2/8  (6) 
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Then  based  on  the  simplified  theory,  the  following  expressions  are  obtained  for  the 
turbulence  intensities  in  the  streamwise  and  crosstream  directions  (Mizukami  et  al.  1992): 

(IT2)  ^/(CuU)  =  (v*2)V2/(CvU)  =  [ed(0/d)2/3/U3]  W  (7) 

where  Cu  =  6.84  ad  Cv  =  4.63.  Results  are  plotted  according  to  Eq.  (7)  in  Fig.  6, 
indicating  excellent  agreement  between  predictions  and  measurements  for  a  1000:1 
variation  of  the  independent  variable  (the  dissipation  factor)  with  both  air  and  water  as 
the  continuous  phase.  Relative  turbulence  intensities  from  this  mechanism  are 
appreciable,  approaching  10%  for  large  values  of  the  dissipation  factor  —  which 
approaches  conditions  typical  of  dense  sprays. 

An  interesting  feature  of  particle-generated  turbulence  is  that  it  exhibits  a  large 
range  of  scales  even  though  the  Reynolds  numbers  of  the  particles  themselves  are 
relatively  low  (<  1000).  This  is  illustrated  in  Fig.  7  by  the  temporal  power  spectral 
densities  of  crosstream  velocity  fluctuations.  These  results  are  for  particles  in  air  but 
results  for  particles  in  liquids  are  similar.  The  measurements  are  accompanied  by  a 
prediction  that  ignores  the  contribution  of  wake  turbulence;  this  limitation  is  not 
significant  because  wake  turbulence  only  makes  a  small  contribution  which  is  at 
frequencies  higher  than  could  be  resolved  by  the  measurements.  Additionally,  sampling 
limitations  prevent  correct  measurements  of  spectra  at  frequencies  higher  than  the  step- 
noise  limit  marked  on  the  plot.  The  unusually  broad  spectrum  is  clearly  evident; 
additionally  the  spectra  decay  according  to  f'1-5  rather  than  f'^/3  which  is  characteristic  of 
the  inertial  range  of  conventional  turbulence.  The  theory  is  in  excellent  agreement  with 
the  measurements  which  helps  explain  these  differences  between  the  spectra;  namely,  the 
turbulence  properties  differ  because  they  result  from  mean  velocity  distributions  in 
randomly-arriving  particle  particles,  with  the  large  range  of  scales  involved  because  flow 
properties  are  strongly  affected  by  wake  mean  velocities. 

Other  spatial  and  temporal  correlations  and  spectra  were  measured  and 
successfully  correlated,  providing  a  rather  complete  description  of  this  novel  form  of 
turbulence,  see  Parthasarathy  and  Faeth  (1990a)  and  Mizukami  et  al.  (1992).  Figure  8 
provides  measured  and  predicted  streamwise  and  crosstream  spatial  integral  scales  as  an 
example.  The  results  are  plotted  as  suggested  by  the  theory  but  reoptimized  to  best  fit  the 
measurements.  The  integral  scales  exhibit  significant  anisotropy,  with  Lux/Luy  =  3.9, 
which  is  typical  of  shear  flows  like  wakes  (Hinze  1975).  The  range  of  scales  also  is 
unusually  large  with  integral  scales  being  150-300  times  Kolmogorov  length  scales, 
which  implies  effective  turbulence  Reynolds  numbers  of  1000-2000  in  terms  of 
conventional  turbulence  properties,  even  though  particle  Reynolds  numbers  generally 
were  ca.  100.  This  behavior  also  follows  because  mean  velocities  of  particle  wakes 
contribute  to  the  turbulence  field. 

3.5  Conclusions 

The  present  investigation  considered  the  properties  of  particle- generated 
turbulence  for  homogeneous  conditions.  The  major  conclusions  are  as  follows: 

1.  Relative  turbulence  intensities  and  integral  scales  could  be  correlated  as  functions  of 
the  dissipation  factor,  edfO/d^/U3.  Relative  turbulence  intensities  from  this 
mechanism  can  be  quite  large,  approaching  10%  for  conditions  comparable  to  dense 
sprays. 

2.  A  number  of  features  of  particle-generated  turbulence  are  similar  to  other 
homogeneous  turbulent  flows:  probability  density  functions  of  velocity  fluctuations 


14 


«d  (0/ d)*/5/u5 


Figure  6.  Streamwise  and  crosstream  r.m.s.  velocity  fluctuations 
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Figure  7.  Temporal  power  spectra  of  crosstream  velocity  fluctuations 
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Figure  8.  Streamwise  and  crosstream  spatial  integral  scales 
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are  Gaussian,  and  normalized  spatial  correlations  and  temporal  spectra  are  relatively 
independent  of  flow  conditions. 

3.  However,  a  number  of  features  of  particle-generated  turbulence  are  distinctly 
different  from  conventional  turbulence  due  to  contributions  from  the  mean 
velocities  of  randomly-arriving  particle  wakes:  the  anisotropy  is  unusually  large 
with  stream  wise  relative  turbulence  intensities  roughly  twice  crosstream  values, 
length  scales  correlate  with  wake  properties  and  are  independent  of  the  average 
spacing  between  particles,  ranges  of  length  and  time  scales  are  unusually  large  even 
though  particle  Reynolds  numbers  are  modest,  and  spectra  decay  at  rates  with 
respect  to  frequency  at  much  lower  rates  than  conventional  turbulence,  e.g.,  f'11 
and  f*1-5  rather  than  f-5#  for  crosstream  and  stream  wise  temporal  spectra. 

4.  A  simplified  model,  based  on  linear  superposition  of  randomly-arriving  particle 
velocity  fields,  was  helpful  for  explaining  many  features  of  the  flow.  However, 
more  information  on  the  character  of  particle  wakes  at  modest  Reynolds  numbers  in 
turbulent  environments  is  needed  to  assess  the  theory  quantitatively  —  work  along 
these  lines  will  be  considered  next. 


4.  SPHERE  WAKES 


4.1  Introduction 

Flow  associated  with  spheres  has  attracted  attention  due  to  numerous  applications, 
e.g.,  dispersed  particle-laden  flows,  sprays  and  rainstorms,  among  others.  The  work  on 
turbulence  generation  by  dispersed  phases  discussed  in  Section  3,  however,  has  shown 
the  need  for  more  information  about  sphere  wakes  at  the  intermediate  sphere  Reynolds 
numbers  (10  <  Re  <  103)  typical  of  drops  in  sprays.  Thus,  the  objective  of  this  phase  of 
the  investigation  was  to  measure  flow  properties  near  spheres  at  intermediate  Reynolds 
numbers,  considering  both  quiescent  and  turbulent  environments. 

A  review  of  past  work  of  the  flow  near  spheres  is  provided  by  Wu  and  Faeth 
(1992).  The  results  of  this  review  indicated  that  there  is  relatively  little  information 
available  about  effects  of  vortex  shedding  on  flow  properties,  and  the  properties  of  wakes 
at  intermediate  Reynolds  numbers,  including  even  whether  these  wakes  are  laminar  or 
turbulent.  Thus,  the  objective  of  the  present  study  was  to  seek  information  to  help  resolve 
these  issues.  Two  experimental  configurations  were  considered,  involving  towed  spheres 
in  a  quiescent  water  bath,  and  spheres  mounted  near  the  axis  of  fully  developed  turbulent 
pipe  flow  in  air.  The  following  description  of  the  research  is  brief,  additional  details  can 
be  found  in  Wu  and  Faeth  (1992a,  b). 

4.2  Experimental  Methods 

Apparatus.  A  sketch  of  the  test  apparatus  for  sphere  flow  in  a  quiescent 
environment  is  illustrated  in  Fig.  9.  The  approach  involved  traversing  test  spheres 
through  a  still  liquid  bath  (415  x  535  x  910  mm)  and  observing  flow  properties  at  the 
center  of  the  bath.  The  sides  and  bottom  of  the  bath  were  insulated  (not  shown  in  the 
figure)  to  minimize  natural  convection  disturbances,  except  for  small  openings  for 
optical  access.  The  bath  liquid  was  water  or  various  glycerol  mixtures  to  provide  a  range 
of  sphere  Reynolds  numbers.  The  test  sphere  was  a  10  mm  diameter  plastic  ball  mounted 
on  taut  125  |xm  diameter  stainless  steel  wire.  The  sphere  was  traversed  at  constant 
velocity  for  a  700  mm  length,  nearly  symmetric  about  the  measuring  location,  using  a 
stepping  motor  driven  linear  positioner.  The  position  of  the  sphere  could  be  traversed 
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Figure  9.  Sketch  ofthe  sphere  wake  apparatus  for  a  quiescent  environment 


19 


normal  to  the  direction  of  the  wire  to  access  different  position  s  n  the  flow,  with  the 
sphere  traverse  itself  provident  access  to  various  streamwise  positions. 

The  test  apparatus  for  spheres  in  a  turbulent  environment  is  illustrated  in  Fig.  10. 
In  this  case,  the  test  sphere  was  mounted  near  the  axis  of  a  fully  developed  turbulent  pipe 
flow  of  air,  within  a  300  mm  diameter  duct.  The  test  sphere  could  be  traversed  in  three 
directions  to  accommodate  rigidly  mounted  instrumentation.  Small  plastic  spheres 
(diameters  of  1.2  -  5.6  mm)  mounted  on  fine  wires  (51  and  127  pm  in  diameter)  in 
tension  were  studied.  Calibration  of  flow  properties  in  the  pipe  indicated  rather  close 
correspondence  of  the  turbulence  properties  to  existing  results  for  fully-developed 
turbulent  pipe  flow,  e.g.  Laufer's  (1954)  results  cited  in  Hinze  (1975)  for  corresponding 
pipe  Reynolds  numbers. 

Instrumentation.  Light  sheet  illuminated  dye  traces  were  used  to  observe  the  flow 
near  the  spheres  in  quiescent  environments.  The  resulting  dye  pattern  was  photographed 
with  exposure  times  of  1  -  4  ms  to  stop  the  motion  of  the  fluid. 

Quantitative  measurements  in  both  flows  were  made  using  laser  velocimetiy  with 
an  arrangement  similar  to  Parthasarathy  and  Faeth  (1990a).  The  results  for  traversing 
spheres  (quiescent  environments)  were  obtained  by  ensemble  averaging  measurements 
from  20  -  120  traverses  at  a  particular  position.  For  spheres  in  turbulent  environments, 
results  were  time-averaged  in  a  conventional  manner. 

4.3  Results  and  Discussion 

Quiescent  Environment.  Flow  visualization  showed  that  the  recirculation  zone  on 
the  downstream  side  of  the  sphere  was  stable  and  symmetric  for  Re  <  200,  stable  and 
unsymmetric  for  200  <  Re  <  280,  and  unstable  with  vortex  shedding  for  Re  >  280.  The 
transition  to  vortex  shedding  was  similar  to  earlier  findings  in  the  literature,  e.g.  transition 
Re  in  the  range  270  -  300. 

Streamwise  mean  velocities  along  the  axis  are  illustrated  in  Fig.  11  for  Re 
extending  up  to  vortex  shedding  conditions.  The  results  show  the  presence  of  three  wake 
regions:  a  fast  decay  wake  region  near  the  sphere  when  vortex  shedding  is  present  (barely 
visible  at  Re  =  280  but  more  pronounced  at  higher  Re),  a  turbulent  wake  region  where  Uc 
~  (x-Xq)'2/3,  and  a  laminar  wake  region  whereuc  ~  (x-Xq)'1  —  the  last  two  regions 
corresponding  to  classical  similarity  behavior  for  turbulent  and  laminar  wakes, 
respectively.  The  fast-decay  wake  region  ends  by  smoothly  transitioning  into  turbulent 
wake  behavior.  Similarly,  the  transition  from  turbulent  to  laminar  wake  behavior  is 
smooth  and  occurs  at  a  wake  Reynolds  number,  Rew,  of  roughly  10.  Thus,  for  low  levels 
of  turbulence  generation,  where  the  continuous  phase  is  nearly  quiescent,  particle  wakes 
largely  exhibit  turbulent  behavior  until  a  final  decay  to  a  laminar  wake. 

Radial  profiles  of  mean  velocities  in  the  important  turbulent  wake  region  are 
plotted  in  terms  of  classical  similarity  parameters,  see  Tennekes  and  Lumley  (1972),  in 
Fig.  12.  Results  of  Uberoi  and  Freymuth  (1970)  for  a  sphere  Reynolds  number.  Re  = 
8600,  and  Chevray  (1968)  for  Re  =  458,000  are  shown  on  the  plots  along  with  present 
results  at  much  lower  Re.  Remarkably,  all  the  results  are  nearly  identical  even  though 
present  measurements  involve  rather  low  wake  Reynolds  numbers. 

The  character  of  the  decay  of  turbulence  near  transition  from  a  turbulent  to  a 
laminar  wake  inillustrated  in  Fig.  13.  Turbulence  intensities  at  the  axis  are  plotted  as  a 
function  of  wake  Reynolds  number  for  present  results  as  well  as  results  from  Uberoi  and 
Freymuth  (1970)  and  Carmodey  (1964)  at  generally  higher  wake  Reynolds  numbers.  For 
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Figure  10.  Sketch  of  the  sphere  wake  apparatus  for  a  turbulent  environment 


Figure  11.  Mean  streamwise  velocities  along  the  wake  axis  for  a  quiescent  environment 
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Figure  12.  Radial  profiles  of  mean  streamwise  velocities  in  the  turbulent 
wake  region  for  a  quiescent  environment 
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Figure  13.  Streamwise  turbulence  intensities  along  the  axis  of  wakes 
in  a  quiescent  environment 
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Rew  >  70,  turbulence  intensities  along  the  axis  are  nearly  constant  with  (uyu)c  in  the 
range  0.85  -  0.92.  At  lower  Rew,  however,  turbulence  intensities  rapidly  decrease  in  a 
final  decay  period.  Based  on  analysis  of  the  final  decay  period  of  axisymmetric  wakes 
(Phillips  1956;  Lee  and  Tan  1967),  behavior  in  this  region  should  be  as  follows  (Wu  and 
Faeth  1992a): 

(u'/u)c  =  C  Rew7/4  (8) 

Present  results  are  in  excellent  agreement  with  Eq.  (8),  with  C  =  1.3  x  10*3  yielding  the 
best  fit  of  the  data.  In  view  of  the  relatively  large  velocity  fluctuations  upon  transition  to 
the  laminar  wake  region  (for  Rew  ~  10)  it  is  remarkable  that  radial  profiles  of  mean 
velocities  in  the  laminar  wake  region  are  inexcellant  agreement  with  classical  laminar 
wake  scaling  (Wu  and  Faeth  1990a).  This  can  be  explained,  however,  by  noting  that 
turbulence  in  the  final  decay  period  is  not  connected  over  the  crossection  of  the  wake; 
instead  it  involves  noninteracting  regions  of  decaying  turbulence  or  turbulence  spots 
(Batchelor  and  Townsend  1948).  Thus,  the  absence  of  connectedness  prevents  mixing  or 
entrainment  by  large-scale  structures  so  that  entrainment  follows  laminar  scaling  with 
mean  velocities  adjusting  to  this  behavior  as  well. 

Turbulent  Environment.  Sphere  wakes  at  moderate  Reynolds  numbers  in  turbulent 
environments  also  exhibited  a  number  of  surprising  features.  The  first  of  these  is 
illustrated  in  Fig.  14,  which  is  a  plot  of  streamwise  velocities  along  the  axis  for  sphere 
Reynolds  numbers  in  the  range  135  -  1560.  In  this  case,  all  the  results  yield  scaling 
satisfying  laminar  wake  similarity  even  though  the  flow  is  clearly  turbulent  (ambient 
turbulence  intensities  are  roughly  4%  and  intensities  within  the  wakes  generally  are 
larger).  This  behavior  can  best  be  explained  by  noting  that  while  the  scaling  is  laminar- 
like,  the  effective  (turbulent)  viscosities  generally  are  2  -50  times  larger  than  the 
molecular  viscosity  over  the  test  range.  Thus,  the  effect  of  the  ambient  turbulence  is  to 
provide  a  turbulent  fluid  whose  effective  viscosity  is  large  enough  so  that  global  scaling 
follows  laminar  scaling  laws. 

The  laminar-like  scaling  of  wakes  in  turbulent  environments  is  clearly  shown  by 
the  radial  profiles  of  mean  streamwise  velocities  plotted  in  Fig.  1 5.  The  scaling  in  both 
the  streamwise  and  crosstream  directions  follows  conventional  laminar  wake  scaling  rules 
(Schlicting  1977)  with  the  effective  turbulent  viscosity  only  varying  with  the  sphere 
Reynolds  number.  As  might  be  anticipated  from  the  results  of  Fig.  14,  laminar  scaling 
provides  an  excellent  representation  of  the  measurements,  well  within  experimental 
uncertainties. 

The  next  important  issue  is  the  behavior  of  the  effective  turbulent  viscosity.  First 
of  all,  it  was  found  that  the  effective  turbulent  viscosity  had  a  single  value  for  a  particular 
sphere  Reynolds  number.  This  behavior  is  illustrated  in  Fig.  16  for  the  present  range  of 
Re  and  sphere  diameters  of  1.2  -  5.6  mm.  The  ratio  vt/v  progressively  increases  with  Re 
but  exhibits  three  ranges:  a  low  Re  regime  for  Re  <  300,  a  transition  regime  for  300  <  Re 
<  600,  and  a  high  Reynolds  number  regime  for  Re  >  600.  The  ratio,  Vj/v  ~  Re  in  the  low 
and  high  Reynolds  number  regimes  with  a  greater  variation  with  Re  in  the  transiton 
regime.  Thus,  the  scaling  in  the  low  and  high  Re  regimes  is  consistent  with  conventional 
eddy  viscosity  ideas  (Tennekes  and  Lumley  1972),  even  though  scaling  of  the  mean 
velocities  themselves  is  laminar-like.  Finally,  the  transition  regime  is  consistent  with  an 
effect  of  onset  of  eddy  shedding.  In  particular,  temporal  spectra  in  the  wake  exhibited  the 
appearance  of  a  spike  when  Re  =  300,  very  similar  to  the  onset  of  eddy  shedding  in 
quiescent  environments.  Furthermore,  the  spike  gradually  disappeared  for  Re  >  600 
suggesting  evolution  into  a  highly  turbulent  regime  consistent  with  the  high  Re  branch  of 
V|/v.  Finally,  it  is  not  surprising  that  the  onset  of  eddy  shedding  is  not  very  different  for 
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Figure  15.  Radial  profiles  of  mean  streamwise  velocities  in  wakes  in 
a  turbulent  environment 
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Figure  16.  Effective  turbulent  \iscosity  for  laminar-like  wakes  in  a 
turbulent  environment 


quiescent  and  turbulent  environments  over  the  present  test  range.  In  particular,  integral 
length  scales  of  the  pipe  flow  are  large  in  comparison  to  sphere  diameters  for  present 
conditions,  ca.  40  mm,  thus,  the  spheres  are  locally  in  a  somewhat  quasisteady 
environment.  Nevertheless,  such  conditions  are  typical  of  spray  environments  where  drop 
dimensions  generally  are  much  smaller  than  integral  scales  (Faeth  1987,  1990). 

4.4  Conclusions 

Aside  from  the  fact  that  the  onset  of  eddy  shedding  always  was  observed  at  Re  ca. 
300,  conclusions  differed  somewhat  for  quiescent  and  turbulent  environments,  as  follows: 


1.  The  wakes  exhibited  fast  decaying  (Re  >  300),  turbulent  and  laminar  wake 
regions,  with  transition  between  the  turbulent  and  laminar  wake  regions  at  Rew  ca 
10. 

2.  Mean  velocities  in  the  turbulent  and  laminar  wake  regions  scaled  according  to 
classical  similarity  results,  even  though  turbulence  intensities  along  the  axis 
varied  in  the  range  10  -  85%,  with  turbulence  spots  extending  into  the  laminar 
region. 

3.  The  main  effect  of  vortex  shedding  was  to  create  a  fast-decay  wake  region  which 
deferred  the  onset  of  the  turbulent  wake  region. 

4.  Turbulence  levels  decayed  according  to  Rew"7/4  as  the  transition  to  a  laminar 
wake  was  approached,  which  is  consistent  with  theoretical  expectations  for  the 
final  decay  period  of  wake  turbulence  (Batchelor  and  Townsend  1948;  Phillips 
1956;  Lee  and  Tan  1967). 


1.  Wakes  in  this  region  were  always  turbulent  but  exhibited  laminar-like  similarity 
scaling  of  mean  velocities  with  effective  turbulent  viscosities  enhancing  mixing 
rates. 

2.  Effective  turbulent  viscosities  were  constant  for  particular  sphere  Reynolds 
numbers  but  progressively  increased  with  increasing  Re,  exhibiting  low  and  high 
Reynolds  number  regimes  separated  by  a  transition  regime  for  300  <  Re  <  600. 

3.  The  transition  regime  was  associated  with  the  onset  of  eddy  shedding  which 
occurred  at  Re  ca.  300,  similar  to  quiescent  environments,  with  the  end  of 
transition  associated  with  the  loss  of  evidence  for  coherent  eddy  shedding  in  the 
temporal  spectra. 

The  present  results  are  limited  because  only  two  ambient  conditions,  quiescent 
and  an  ambient  turbulence  intensity  of  ca.  4%,  were  considered.  Thus,  behavior  at  other 
conditions  must  be  known  before  the  results  can  be  applied  to  analysis  of  turbulence 
generation.  Information  concerning  scalar  mixing  in  sphere  wakes  at  modest  Reynolds 
numbers  also  is  required  to  understand  the  scalar  turbulent  field  of  dense  sprays.  Both 
issues  are  being  addressed  during  current  work  in  this  laboratory. 
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Figure  2.  Drop  breakup  Umes  as  a  function  of  We  and  Oh. 


5?  *  « 
“  S  o  ^ 


I  2  i 

5  -o  ^  i 

§  5  1 1  i 

*i  |«] 

Sm  £5 

e  a  -o  “ 

„  ■*>  5  o 
iS=>.J 

S  c  “4  i 

1  I !  |  i 

55  i?  S  **  < 

■g  a  H J. 

..  g  a.  S  ■ 

i**!; 

E-«SS' 
g  3|  2, 

2  S  I  ?  1 

2  *  3  | 

ax«  3 
E"a«; 

3  -g  S  S 
«  S-o  « 

ft*  g 

»  =  o  <2 

=  K  a  § 
«  «■?  2 
■-  -p  «  2 

555  I 
~  %  Z.  jj 

1 !  S  4 

Sica, 


|!|  | 


a«s  3  s 

s^is  a 


-  -6  o  a 
■J  S  -5  g 
2  l!  a 
i  3.  o  -5 


■wSlfi  al 

1  I  gA*! * 
o  c-  if  - 1  ° 

C  O  u  2  E  -2  r 

2  £  5  “  c  £  Ji 
rt  o  :r  J2  J  rj  *= 

iiiiiii 

m  |  “z  S  o  s 

|  S<*-~  -5 

vr:|vi* 

ihh^i 

Lo^vJi 

«  «  I  2-4-3  ° 

|a J  § :££ 

£  -g  2  v  -2  3“  a 

|  =  *  £  3  -2  * 

SgsUgB 

5  £  ■'»  a  o  2  | 

s  3  E  B  *■  E  “ 

■>3  s  5  2i-s<i 
,:Sh3J!;ssI 
j  5.  .  2  =  e  -  5  V 
j  3  §-s  *  Ij  1 3 

1 1  3  1 13!  §■! 

5  I  -f  S  s  «  S  3  * 

i  I  °  v-55  s-6! 

'EsIibssS 

;|lf  U\U 

^illf  §|  I 

i  3  c  o  Cl  -5  > 


§  o  s  s' *  o  a 
s  ^  i  Ejs  ?J 
«s  s  s i 
S*!  H  »* 

|!|325 • 

(9  2  "Z  *  £  9“ 
§  3 i  f  V 

•X  -  ^  **  V  '■*  «J 

S3 3  5  3  §* 

| ." !  .s !  i| 

S  ~  t:  a  &  si  2 

O  5=8  J. 

S l! !* 

£  3  3  2P  3  — 

■8  1!  3  ■*  £  3  * 

Cl  -C  ^  L3  ^  «  S* 

3s|«j?=. 
~  2  s  5  s  “  = 
a  =  |  2  “  'S  15  ■ 

2,-2  5  if  .j  >■ 

;  £  g  5  s*  sr  -3  £  . 

i  °  2  s^j  !  °- 

iimiii. 

ililtilli 

2  .  o.  a  -=  -f  c 


Ilf 

*  g : 


43 


Ill 

C  an  E. 


I  >“8  O  ?  -g 

Sieges 

*  o°  8  I  1 1 

111  3“  8 

51  i-* »« 

p£  «  s  ‘  i 
a  «  3  1  s  “ 

*!§fif 

§  g  £L  -<?  I  2 
;f^ss 
l  ~  g  o  ^  ■£ 

I  r*  £  ,/>  r-  *• 

a  E  —  o 
Sou-' 

Ulifl 

!«c»2" 
c  e  a.<£:  ^ 

S  S  .2  o  c  S 

3  ^  I  *  %  a 

j 2  if«i 

IllSfc 

J  X  u  5  2 

• :  I  s  *  § 

I s  *  I !  J 

lUUl 

i * ®  ^  1 " 

*  £  c  -  5»  o 

*  -  o  -  —  .3 

*  «-  s  -c  !a  n 
\  °  al’S 

.  *  3  fi  «*  * 

I  2  w  ®*  w  S 


1  'I  { p , 

is=li 

«•  *2  w  -fi  2 

£  3r  ag 
.S  iC  “  ^  i« 
x  U  u.  -a 

2  s  3  .£• 

S  f-4  t  Q  « 

$■  n  s  t  = 

‘  j  ,5  « 
,o3  o.1^  « 
'uvSs: 
c  -  -a  a  s 

—  a  u  * 
§•8  “  2  .§ 

S  -5  o  a  -5 
g  £  c  «•  O 

O  y>  v 

■tf  -  c  §  » 

*«2  5  5  » 
u  <a  j:  x 
o  .r  73  E  a. 
w  a  >>  o 
u-  C  <2  ^ 

.3^*^ 1 

25  o  £  o  a 

3  a.  E  £  £ 

J8--8.il 

^  Ss?| 

£  £  C  8  v  _4 

v*i«;5 

d  S  «S  5  " 
S  6  3  * 

silsi! 

i  ^  ^  y  JJ  ti 
“  5  ■§  =  ^  c 

g  S  4  *  3  S 

S  -S.J  s  s  s 

Is  II 2  ? 


>  s,  a  s 

'1  “■* 

•  r  ii  M 

5.2-S.2 

>  «  Js  S» 

»  oo<£  •£ 

I  §  as 

:  «  S  -8 

)  c  c 

ills. 

1 E  2 


2  o  v  J!  •>  5 

gig 

6  5  s  s  8  o. 


.O  D.  —  O  O 
w.  O  ^  C  C 

ii  i  * » 

■sf  111 
ill  °  § 
1 11 !  j 

>.  3  §  ]^ 

S  -g  2  2 

|-  §  s-g 

v-as-C" 

Od>  C<«  s 

»i  '  5  o  .2 
—  -5  c  <3 
y  J  -  O  2 
y  3s  2  S  v 

S-  i  3  2 

J*  §|  * 

;  *  1  - 1 

|*i!i 

^  ^  «  r»i  a 

a  *  1 S  § 
S  5“  “ 
2  o£ja-= 

»-.;i  s 
I’siij 

s 

:  *  *o  s'  T  c 
:  «  3  2  M  o 

1  c-  £5  c 
J-  0  g  *®  vi  <2 
»  *2  Ji  c  o  c 
\  |  5  §  5  .2 
:  s  1 1 5  i 
1**131 

*  A  -a  i  "o 


44 


The  present  results  yielded  btmodal  behavior  for  shear  breakup  (figure  8)  where  departure  from 
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CUMULATIVE  VOLUME  PERCENTAGE 

Figure  8  Distribution  of  drop  diameters  after  shear  breakup  Figure  **  Sketch  of  the  shear  breakup  process 
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Appendix  C:  Parthesarathy  and  Faeth  (1990a) 
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Turbulence  modulation  in  homogeneous  dilute 
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range  would  be  helpful  in  resolving  this  issue.  LV  burst  processor  introduced  step-noise'  caused  by  the  sample-and-hold  signals  of 

The  p  d  f  a  of  streamwise  and  cross-stream  velocity  fluctuations  are  plotted  as  a  the  processor.  Adrian  &  Yao  (1987)  have  shown  that  effects  of  step-noise  are 

function  of  normalized  variables  in  figures  ti  and  7.  Within  experimental  observed  at  frequencies  roughly  one-tenth  of  the  frequency  corresponding  to  the 

uncertainties,  measurements  are  generally  independent  of  |>artiele  size  and  loading.  mean  data  rate.  The  characteristics  of  step-noise  are  the  appearance  of  a  noise  band 

Predictions  based  on  a  Gaussian  p.d.f.  which  is  the  form  expected  from  the  theory.  having  a  constant  spectral  amplitude  followed  by  a  second -order  low- pass  filter  effect 
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this  condition  provides  a  reasonable  indication  of  effects  of  turbulent  dispersion.  Ill  dis|iersjon  for  the  0.5  nun  diameter  particles  In  fact.  iiotntK  that  the  ratio  of 

eontrast.  eross  stream  velocity  fluctuations  due  to  self-induced  particle  motion  are  (e"’)l/(r',|l  is  roughly  eonstant.  indicating  similar  particle  res|snise  to  cross  stream 

generally  larger  than  those  observed  in  the  bath  for  the  1.0  and  2.0mm  diameter  liquid  velocity  fluctuations  over  the  |iresent  test  range  (see  table  II)  this  iuerease 

particles  so  that  these  test  conditions  are  of  questionable  value  for  gaining  simply  follows  the  <  •  variation  of  liquid  velocity  Huetuatiniis  Surprisingly  the  cross- 

information  about  turbulent  dis|*-rsinn  A  curious  phenomena  is  that  cross-stream  stream  particle  velocity  Hurt  nations  also  exhibit  significant  increases  with  increasing 
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Appendix  E:  Mizukami  (1992) 
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Appendix  F:  Wu  and  Faeth  (1992) 
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Graduate  Assistant,  Department  of  Aerospace  Engineering. 

*  Professor,  Department  of  Aerospace  Engineering.  Fellow  ALAA. 
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periphery  of  the  sphere,  giving  the  appearance  of  a  vortex  street  passing  into  the  sphere  u,lu/„n  ..  ,  .  ,  , 

between  these  regimes  and  the  scaling  of  turbulence  properties  at  intermediate  wake 

wake  when  viewed  in  crossection.15'16  These  results  were  obtained  from  flow 
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traverse*  at  a  particular  radial  position  and  sphere  Reynolds  number.  The  number  of  kinematic  viscosities  were  compensated  by  adjusting  the  sphere  velocity  to  obtain  the 

traverses  was  selected  to  control  the  experimental  uncertainties  and  were  large  in  regions  desired  sphere  Reynolds  number. 


Drag  coefficients  for  the  various  sphere  Reynolds  numbers  generally  were  The  similarity  results  of  Eqs.(2M4)  yield  Ue-Kx-Xo)'^3,  and  Re*,  ~  (x- 

compuied  from  correlations  (hie  to  Putnam28  and  KUrten  et  al.29  The  one  exception  was  Re  *»  fully-devetoped  turbulent  wakes  far  from  the  source  —  the  last  implying  eventual 

*  280,  which  was  at  the  onset  of  vortex  shedding  where  existing  drag  correlations  were  transition  to  a  laminar  wake  as  Rew  becomes  small  if  the  surroundings  are  quiescent.21 

suspect;  therefore,  Co  at  this  condition  was  found  from  the  measured  velocity  defect  in  the 
wake.  Values  of  the  initial  momentum  thickness  of  the  wakes  then  were  computed  as  1  -aminar  w»in- 


velocity  and  characteristic  wake  width  then  becomes:  growth  of  the  width  of  the  wake,  are  larger  for  laminar  than  turbulent  wakes. 
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be  discussed  later  in  more  detail,  based  on  visualizations  extending  farther  into  the  wake.  shedding  process,16  even  at  the  highest  sphere  Reynolds  number  illustrated  in  Fig.  4. 

Additionally,  the  vortex-like  (swirling)  features  seen  near  the  outer  boundary  of  the  dye 
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considered.  The  parameter  (U0)i/Us  »  2.23  for  exact  agreement  with  Eq.  (2):  present  a  maximum  in  the  screamwise  distance  required  to  develop  a  turbulent  wake  structure  for 

measurements  agree  with  this  estimate  for  Re  £  400,  with  an  average  value  and  standard  Re  <  103. 
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The  radial  profiles  of  mean  streamwise  velocities  in  the  turbulent  wake  region  are  28°-4000  (present  measuremen“  were  not  feasible  at  lower  values  of  Rc  because  excessive 

illustrated  in  Fig.  8.  In  addition  to  present  measurements,  results  from  Uberoi  and  numbers  of  mnverscs  were  ,0  establlsh  reliab,e  values  ofC‘  due  ,0  low  si*na,t<>- 

Freymuth19  for  a  sphere  having  a  Re  =  8600  and  Chevray20  for  a  siender  body  having  a  Re  noise  ratios)'  Results  for  lhc  fuU  nn*e  of  *»  for  the  measurements  « illustra,cd  wuh  the 

portion  in  the  turbulent  wake  region  (or  the  fully-developed  turbulent  wake  region 
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constant  in  the  self-preserving  wake  region  (this  latter  portion  corresponding  to  behavior 
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Fig.  1 1  Correlation  of  radial  profiles  of  mean  streamwise  velocities  in  the  laminar 
wake  region  (Re  =  90-280). 
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